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6.  A  One  Year  Technical  Research  Summary: 


-“Error-Trellis  Syndrome  Decoding  for  Convolutional  Codes." 


A  new  error-trellis  syndrome  decoding  scheme  for  CCs  is  developed.  It 
is  demonstrated  that  the  real  advantage  of  error-trellis  decoding  over 
both  Viterbi  and  sequential  decoding  of  CCs  is  the  reduction  of  the 
number  of  states  and  transitions  between  any  two  frames. 


"CSI  Architecture  for  Algebraic  Syndrome  Decoding  of  Dual-K 
Convolutional  Codes.’ 


An  algebraic  syndrome  decoder  is  developed  to  find  the  best  estimated 
message  sequence  for  dual-K  CCs  without  finding  minimum-error  paths 
in  an  error-trellis  diagram.  The  advantage  of  this  algebraic  syndrome 
decoder  over  an  error-trellis  decoder  of  the  dual-K  CCs  is  that  the 
message  sequence  can  be  corrected  without  a  necessity  for  storing  a 
large  number  of  states  or  paths  in  a  constraint  length  of  the  error  trellis 
diagram.  Finally,  a  LSI  chip  is  developed  to  realize  this  algorithm. 


-"VLSI  Design  of  a  Pipeline  Algebraic  Syndrome  Decoder." 


A  new  VLSI  architecture  is  developed  for  the  Algebraic  Syndrome 
decoder.  The  advantage  of  this  new  architecture  is  that  a  substantial 
reduction  in  the  number  of  transistors  is  accomplished. 


-"Searching  High-Rate  Systematic  Optimum  Distance  Convolutional 
Codes." 


Some  high-rate  systematic  optimum  distance  convolutional  codes  are 
being  searched  with  rates  up  to  7/8  and  of  constraint  length  up  to  15. 
These  codes  can  be  efficiently  decoded  using  error-trellis  syndrome 
decoding. 
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-'Decoder  Performance  Simulation.' 

Decoder  performance  simulation  is  accomplished  both  for  error-trellis 
decoding  of  convolutional  codes  and  for  algebraic  syndrome  decoding  of 
dual-K  convolutional  codes. 

7.  Published  Papers  and  Abstracts. 

-I.  S.  Reed  and  T.  K.  Truong,  "Sequential  Syndrome  Decoding  Techniques 
for  Convolutional  Codes,"  submitted  to  IEE  Proceedings,  pt.  E. 

Abstract  -  This  paper  reviews  previous  studies  (Refs.  1  and 
2)  of  the  algebraic  structure  of  convolutional  codes  and 
extends  those  studies  to  apply  to  sequential  syndrome 
decoding.  These  concepts  are  then  used  to  realize  by  example 
actual  sequential  decoding,  using  the  stack  algorithm. 

-J.  F.  V-ing,  I.  S.  Reed,  T.  K.  Truong,  J.  Sun  and  J.  Y.  Lee,  "LSI  Architecture 
for  Algebraic  Syndrome  Decoding  of  Dual-K  Convolutional  Codes," 
submitted  to  IEE  Proceedings,  pt.  E. 

Abstract  -  In  this  paper,  algebraic  syndrome  decoders  are 
developed  which  extend  the  early  syndrome  decoders  of 
certain  convolutional  codes  such  as  the  Wyner-Ash  code. 
Specifically  syndrome  decoders  are  designed  to  decode  both 
the  rate  1/2  and  1/3,  dual-k,  nonsystematic  convolutional 
codes  (CCs).  Also  the  LSI  architectures  of  these  decoders  are 
presented.  Further,  it  is  demonstrated  that  such  decoders  can 
be  realized  readily  on  a  single  chip  with  CMOS  technology. 

The  advantage  of  this  algebraic  syndrome  decoder  over 
error-trellis  decoding  of  dual-k  CCs  is  that  the  message 
sequence  can  be  corrected  without  the  necessity  for  storing  a 
number  of  states  or  paths  in  a  constraint  length  of  the  error 
trellis  diagram. 

-J.  M.  Jensen  and  I.  S.  Reed.  "Bounded  Distance  Coset  Decoding  of 
Convolutional  Codes,"  IEE  Proceedings,  vol.  133,  pt.  F,  no.  5,  August  1986. 

Abstract  -  This  paper  presents  a  maximum  likelihood 
consistent  bounded  distance  decoding  algorithm  for 
convolutional  codes.  The  algorithm  correctly  decodes  all  error 
sequences  which  fall  within  the  error  correcting  sphere.  A 
class  of  codes  is  defined,  in  which  the  decoder  exploits  the 
fact  that  only  certain  error  sequences  need  to  be  corrected. 

For  these  codes  the  decoding  is  based  on  a  reduced  encoder 
state  diagram.  Thus  only  a  subset  of  the  trellis  or  tree  has  to 
be  searched  in  order  to  find  the  error  pattern.  An  exact 
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characterization  of  the  reduced  state  diagram  is  given  in  this 
paper  along  with  an  example. 

-H.  M.  Shao,  T.  K.  Truong,  I.  S.  Hsu,  L.  J.  Oeutsch  and  I.  S.  Reed,  "A  Single 
Chip  VLSI  Reed-Solomon  Decoder,*  Proc.  Int'l.  Conf.  on  Acoustics, 
Speech  and  Signal  Processing,  Tokyo,  Japan,  April  7-11,  1986. 

Abstract  -  A  new  VLSI  design  of  a  pipeline  Reed-Solomon 
decoder  is  presented.  The  transform  decoding  technique  used 
in  a  previous  design  is  replaced  by  a  simple  time  domain 
algorithm.  A  new  architecture  which  realizes  such  algorithm 
permits  efficient  pipeline  processing  with  a  minimum  of 
circuits.  A  systolic  array  is  also  developed  to  perform  erasure 
corrections  in  the  new  design.  A  modified  form  of  Euclid's 
algorithm  is  developed  with  a  new  architecture  which 
maintains  a  real-time  throughput  rate  with  less  transistors. 

Such  improvements  result  in  both  an  enhanced  capability  and 
significant  reduction  in  silicon  area,  thereby  making  it  possible 
to  build  a  pipeline  (255,223)  RS  decoder  on  a  single  VLSI  chip. 

-J.  F.  Wang,  I.  S.  Reed,  T.  K.  Truong  and  J.  Sun,  "Algebraic  Syndrome 
Decoding  of  Dual-K  Convolutional  Codes,"  to  be  submitted  for 
publication  soon. 

In  this  paper,  algebraic  syndrome  decoders  are  developed 
which  extend  the  early  syndrome  decoders  of  high  rate 
convolutional  codes  such  as  the  Wyner-Ash  code.  In  this 
paper,  syndrome  decoders  are  designed  to  decode  the  rate  1/n 
dual-k  nonsystematic  convolutional  codes.  The  advantage  of 
the  algebraic  syndrome  decoders  over  error-trellis  decoding  of 
dual-k  convolutional  codes  is  that  the  message  sequence  can 
be  corrected  without  the  necessity  of  storing  a  large  number 
of  states  or  paths  in  a  constraint  length  of  the  error  trellis 
diagrams. 

-I.  S.  Reed,  I.  S.  Hsu,  J.  M.  Jensen  and  T.  K.  Truong,  "The  VLSI  Design  of 
an  Error-Trellis  Syndrome  Decoding  for  Certain  Convolutional  Codes," 
IEEE  Trans,  on  Computers,  vol.  C-35,  no.  9,  pp.  781-789,  September  1986. 

A  recursive  algorithm  using  the  error-trellis  decoding 
technique  is  developed  to  decode  certain  convolutional  codes, 
such  as  dual-k  convolutional  code.  It  is  demonstrated  that 
such  a  decoder  can  be  realized  readily  on  a  single  chip  with 
NMOS  technology. 
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ABSTRACT 

In  this  report,  algebraic  syndrome  decoders  are  developed  which  extend  the  early  syndrome 
decoders  of  certain  convolutional  codes  such  as  the  Wyner-Ash  code.  Specifically  syndrome  decoders 
are  designed  to  decode  both  the  rate  1/2  and  1/3,  dual-k,  nonsystematic  convolutional  codes  (CCs). 
Also  the  LSI  architectures  of  these  decoders  are  presented.  Further  it  is  demonstrated  that  such 
decoders  can  be  realized  readily  on  a  single  chip  with  CMOS  technology. 

The  advantage  of  this  algebraic  syndrome  decoder  over  error-trellis  decoding  of  dual-k  CCs  is 
that  the  message  sequence  can  be  corrected  without  the  necessity  for  storing  a  number  of  states  or 
paths  in  a  constraint  length  of  the  error  trellis  diagram. 


I.  INTRODUCTION 


Recently  the  authors  [1,2,3]  developed  a  new  erTor-trellis  syndrome  decoding  scheme  for  CCs. 
This  method  involved  finding  minimum-error  paths  in  an  error-trellis  diagram,  it  was  shown  [1,2,3] 
that  the  computation  of  the  error  trellis  could  be  accomplished  by  finding  the  solution  of  a  syndrome 
equation  explicitly  in  terms  of  the  actual  error  sequences.  The  real  advantage  of  error-trellis  decoding 
over  Viterbi  coding-trellis  decoding  of  CCs  is  the  reduction  of  the  number  of  states  and  transitions 
between  any  two  frames. 

It  is  shown  (e.g.,  see  [4])  that  a  simple  logic  design  for  a  syndrome  decoder  can  be  found  for  a 
rate  3/4,  one-error-correcting  systematic  Wyner-Ash  code.  In  this  paper,  it  is  shown  that  these  early 
syndrome  decoders  such  as  that  used  for  the  Wyner-Ash  code  can  be  extended  to  decode  rate  1  In, 
dual-k  CCs  for  n  =  2,3.  It  is  well  known,  see  [5],  that  a  rate  1  In  dual-k  CCs  is  capable  of  correcting 
only  r  errors  in  two  codeword  frames  where  t  =  n-1.  From  this  fact,  it  is  shown  here  that  an  alge¬ 
braic  syndrome  decoder  can  be  developed  to  find  the  best  estimated  message  sequence  v  for  dual-k 
CCs  without  finding  minimum-error  paths  in  an  error-trellis  diagram.  In  other  words,  if  no  more  than 
t  symbol  errors  occur  in  two  codeword  frames,  an  algebraic  syndrome  decoder  can  be  developed 
which  recursively  determines  the  best  next  corrected  message  symbol  v,  from  only  the  current  best 
estimated  symbol  v(_!  and  the  received  sequence  at  time  frames  i  and  i+l. 

Next,  an  LSI  architecture  is  developed  to  realize  these  new  algebraic  syndrome  decoders  for  both 
rate  1/2  and  1/2,  dual-3  CCs.  The  designs  of  these  decoders  are  regular,  simple  and  therefore  naturally 
suitable  for  LSI  implementation. 

n.  PROPERTIES  OF  CONVOLUTIONAL  CODE 

In  order  to  systematically  develop  an  algebraic  syndrome  decoder,  certain  properties  of  a  convo¬ 


lutional  code  are  needed. 


First  let  the  information  of  message  sequence,  the  input  to  the  CC,  be  represented  by 

x(D)  =  [xx(DU2(D)  xk(D)],  (la) 

where 

*,(£>)=  I  l<j<>k  (lb) 

i-0 

are  elements  in  F[D],  the  ring  of  polynomials  in  the  unit  delay  operator  D  over  F  =  GF{q),  a  Galois 
field,  with  q ,  a  power  of  a  prime  integer.  Vector  x(D )  is  a  generating  function  in  D  of  the  input  mes¬ 
sage  sequence  x  =  [xq,  •  •  x;,  ■  ■  •  ],  where  x;  =  [x;1,  .  .  .  ,xjk\  is  a  vector  belonging  to  VF(k),  the  k- 
dimensional  vector  space  over  F .  x(D)  is  sometimes  called  a  D-transform  of  the  message  or  informa¬ 
tion  sequence  x.  The  k  component  vector  xj  in  x  is  called  the  information  frame  at  stage  or  frame 
time  j . 

In  a  similar  manner,  the  output  sequence  is 

y(D)  =  [y,(D) . y„(D)],  (2) 

where  y,(D)  e  F[D],  15/  <  n.  Vector  y(D)  is  the  D-transform  of  output  coded  sequence 
1  =  (Zo<  ■  •  where  =  [y7i,  •  •  •  ,yjn]  belongs  to  Vn{F).  The  n-vector  y}  is  called  the  jlh 

codeword  frame  of  code  sequence  y . 

The  information  and  code  sequence  of  an  (n,k)  convolutional  code  are  linearly  related  by  a 
k  x  n ,  rank  k ,  generator  matrix  G  (D )  of  polynomial  elements  in  F  [ D  ],  as  follows: 

y(D)  =  x(D)G{D).  (3) 

The  maximum  degree  m  of  the  polynomial  elements  of  G  {D )  in  D  is  called  the  memory,  and  the  con¬ 
straint  length  L  is  defined  as  L  =  m+1. 

The  free  distance  of  a  CC  is  defined  by 
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dfne  ~  Min  Wh(z(D)),  (4) 

y(D)*Q 

where  WH(£(D))  is  the  cumulative  Hamming  weight  of  the  coefficients  of  DJ  for  all  j  >  0,  where 
Vj  is  the  j,h  codeword  frame.  Note  that  the  computation  of  dtne  requires  at  least  L  codeword  frames 
for  all  codes  of  practical  interest. 

To  avoid  catastrophic  error  propagation,  G  (D )  is  assumed  to  be  in  a  format  of  a  basic  encoder 
[6],  The  Smith  normal  form  of  a  basic  encoder  [2]  is 


G{D)  ~  A{D  )[/*,O]B(0). 


(5) 


where  A  ( D )  and  B  (D )  are,  respectively,  k  x  k  and  n  x  n  invertible  matrices  over  F  [£>  ]  and  Ik  is  a 
k  x  k  identify  matrix. 

In  Eq.  (5),  let  matrix  B  [D  ]  be  partitioned  as 


B(D)  = 


B\(D) 

B2{D) 


where  B\(D)  consists  of  the  first  k  rows  of  B(D)  and  B2(D)  consists  of  the  last  n-k  rows  of  B(D). 
Similarly,  let 

fl(D)-1  =  (S,(D)^2(D)], 

where  B  ,(D )  consists  of  the  first  k  columns  of  B  (D  )-1  and  B  2(D )  consists  of  the  last  n  -k  columns  of 
B(D)~l.  Since  B  (D)  B (D)~{  =  /„,  the  following  identities  evidently  hold: 

B \{D  )  B )  =  /*,  B]{D)  B2(D)  =  0 

(6) 

52(D)5,(D)  =  0,  B2(D)B2(D)  =  l„_k 

A  parity-check  matrix  H(D )  is  defined  to  be  an  (n-k)  x  n  matrix  of  rank  ( n-k  ),  satisfying 


G(D)Ht(D)  =  0, 


(7) 


From  Eqs.  (5),  (6),  and  (7),  it  is  seen  next  that 


H(D)  =  BiT(D) 


has  the  properties  of  a  parity-check  matrix  H(D )  associated  with  G  ( D ). 


By  Eq.  (3),  the  CC  generated  by  G(D)  is  the  set 


C  =  (2(D)  =  [y,(D) . y„(D)]|y(D)  =  £(D)G(D)}. 


It  is  now  shown  also  that 

C  =  {1(D)  =  [yx(D) . y„(D)]|y(D)//r(D)  =  0},  (10) 

where  //(D)  is  given  in  Eq.  (8).  To  see  this,  denote  the  right  side  of  Eq.  (10)  by  C Clearly  an  ele¬ 
ment  of  C ,  as  given  in  Eq.  (9),  belongs  to  Cw  and  hence  C  £  CH. 

Next  suppose  y_\{D)  is  an  element  of  CH,  i.e.  by  Eqs.  (8)  and  (10), 

lx(D)HT{D)  =  ll(D)B2(D)  =  0. 

But,  by  definition,  B2(D)  consists  of  the  last  ( n-k )  columns  of  B~l{D),  so  that 


B2(D)  =  B~\D)  , 

1  n-k 


where  ”0"  denotes  a  block  of  k  rows  of  zeroes  and  /„_*  is  the  {n-k)  row  identity  matrix.  Thus. 
V](D)  satisfies  the  equation 


lx(D)B~\D)  0  =0. 

1  n-k 


The  most  general  solution  of  this  equation  for  y  j(D  )B~\D )  is 

lx(D)B-\D)  =  (X,(D ) . xt(D),0  0]  =  [X(D),0] 

where  x;(D)  for  1  <  ;  <  k  can  be  chosen  to  be  any  arbitrary  elements  of  F[D}.  Solving  for  yi(D) 


yields,  finally,  by  Eq.  (5), 


y,(D)  =  T (D  )[/*,03fl(D )  =  X(D  )A~'(D)G  (D ), 
which  belongs  to  C ,  as  given  in  Eq.  (9).  Thus  Q  cC  and  Eq.  (10)  is  proved. 

The  fact  that  the  CCs  given  by  the  set  C  in  Eq.  (9)  can  be  characterized  by  Eq.  (10)  is  used  in 
the  following  section  to  find  the  coset  of  solutions  to  the  syndrome  equation. 

HI.  METHOD  OF  ALGEBRAIC  SYNDROME  DECODING 

Let  y(£> )  in  Eq.  (3)  be  transmitted  and  z(D )  be  received.  Then, 

z(D)  =  1(D)  +  e(D),  (12) 

where  e(D)  is  the  D-transform  of  the  error  sequence.  By  Eqs.  (12)  and  (7),  the  syndrome  of  the 
received  sequence  is 

1(D)  =  z(D)HnD)  =  [i(D)  +  e(D)]HT(D) 

=  e(D)HT(D),  (13) 

or  its  equivalent, 

(e(D)  -  z(D))HT (D)  =  0,  (14) 

for  all  solutions  e(D). 

By  Eqs.  (10)  and  (9),  the  term  (e(D )  -  z(D ))  in  Eq.  (14)  must  be  some  code  sequence 
v  (D  )G  (D ).  Hence,  the  most  general  solution  of  the  syndrome  equation,  Eq.  (14),  is 

e(D)  =  z(D )  +  v(D)G(D ).  (15) 

where  v(D)  is  the  D-transform  of  an  arbitrary  message-like  sequence  y  =  [y0,  .  .  ,\>j,  ■  ]  of  k- 

vectors  v;  e  vt(F). 


Equation  (15)  shows  that  the  most  general  solution  of  the  syndrome  equation,  Eq.  (13),  for  e(D) 
is  the  coset 

Ct  =  (e(D )  =  z(D)  +  v(D)G(D)ly(D)  =  [v  ,(D ) . v*(D )]} 

of  code  C,  defined  by  either  Eq.  (9)  or  Eq.  (10).  A  minimization  of  the  Hamming  weights  over  all 
elements  of  coset  C2  yields  the  standard  minimum-error  solution  for  message  y(D).  Efficient  methods 
for  achieving  this  minimization  include  the  Viterbi  algorithm  and  all  sequential  decoding  methods  for 
convolutional  codes. 

The  difficulty  with  the  standard  decoding  methods  of  CCs,  i.e.,  Viterbi  or  sequential  decoding,  is 
the  need  to  consider  a  sometimes  prohibitively  large  number  of  states  and  paths  in  the  decoding  trellis. 
Such  minimum-weight,  path-finding  decoding  method  does  not  take  the  advantage  of  the  limited  error- 
correcting  capability  which  one  might  expect  could  reduce  complexity  of  the  decoder.  It  is  shown  in 
the  next  section  that  such  reduced  complexity  is  possible  by  using  Eq.  (15)  to  systematically  develop 
algebraic  syndrome  decoders. 


IV.  ALGEBRAIC  DECODING  OF  DUAL-K  CCs 

Dual-k  convolutional  codes  are  of  rate  1  In  of  memory  m  -  1,  and  with  symbols  in  the  finite  or 
Galois  field  GF( 2*)  (see  Odenwalder’s  paper  [5]).  The  generating  matrix  G  has  the  following  form, 
namely. 


G  = 


G0  Ci 
G0  G, 


(16a) 


where  G0=  [1,1,1.  ,1]  and  Gx  =  fg;i,gi2.  •  .gi„l  with  gij  *  0  and  gtj  €  GF( 2*)  and  the  gtj' s 

are  all  distinct,  for  1  <  j  5  n .  Thus  by  (3) 


G(D)  =  G0  +  G,D 


(16b) 


where  D  is  the  unit  delay  operator. 

From  the  above  definition  of  a  dual-k  CC,  it  is  known  in  [5]  or  can  be  verified  that  the  minimum 
distance  of  the  code  is  =  (2n-l)  and  the  free  distance  is  dfree  =  In.  Hence  if  no  more  than  t 
symbol  errors  occur  in  the  first  2  codeword  frames  and  2r+l  <,  n  =  2n-l  or  t  <  n-1,  then  those 
errors  which  occur  in  the  first  frame  can  be  corrected.  In  other  words,  the  dual-k  CC  is  a  t-error-per- 
blocklength-correcting  CC,  where  t  =  [(<d min- 1  )^2]  and  [x]  denotes  the  greatest  integer  less  than  x. 
Note  that  blocklength  is  equal  to  n(m+l)  =  nL  symbols  where  m  is  the  memory  and  L  is  the  con¬ 
straint  length. 

The  following  sections  present  the  algebraic  syndrome  decoding  algorithm  for  the  rate  1/2  and 
1/3  dual-k  convolutional  codes. 

V.  ALGEBRAIC  SYNDROME  DECODING  OF  A  RATE  1/2,  DUAL-K  CC 

Let  the  Galois  field  GF( 23)  be  generated  by  the  3rd  degree  irreducible  polynomial 
p(x)  =  x3+x2+l,  over  GF ( 2).  If  a  is  a  root  of  p(x),  then  a,  a2,  a3  =  1+a2,  a4  =  1+a+a2,  a5  =  1-KX, 
a6  =  a+a2,  a7  =  1  and  0  are  the  eight  elements  of  GF( 23).  The  generating  matrix  of  type  (16b)  for  a 
rate  1/2,  dual-3  CC  is  given  by 

G  (D  )  =  [l,a(]  +  [ I ,a 2 J D  =  [  I+fl  \D , I+a J. 

where  <3[  *■  a2t  a2  *  0.  In  Eq.  (15),  let 

y(0)=£v,D'  (17) 

1*0 

£(fl)=  £  leu'‘2,V>‘ 

and 


%  V  \  \ 


(18) 


l(D)=  £  [r u,z2i]D 
1=0 


A  substitution  of  Eqs.  (17),  (18),  (19)  into  Eq.  (15)  yields 


Z  [e1(,e2l]D‘  =  £  v,D‘G(D)+  £  [zu,z2i]Dl 

t  =0  i=*0  I  =0 


Z  =  Z  lv«  +  alv»-l-  vi  +  *2vi-llC>‘  +  Z  [z!oz2JP‘ 

1=0  i=0  1=0 


Thus, 


v,  +a,v,_[  =  z,,-  +eu  (21 

vi  +  «2VM  =  z2i  +  «2.  (22 

Solve  v,  and  v,_j  in  terms  of  zh,  22,  ,  e1(  and  e2l,  one  obtains, 

v,-i  =  (a  1  +  «2)“1K*ii  +  z2i)  +  («1.  +  «2i)]  (23 

V,  =  (ai  +  a2r1[(a2zi,  +  <*iz2i)  +  ( a2eu  +  a,e2l)]  (24 

where  v,_!  and  v,  constitute  the  current  message  symbol  and  next  message  symbol  in  the  input  infor 

mation  sequence. 

In  the  above  notation,  the  following  results  can  be  established. 


Result  1: 


Let  the  error  correcting  capacity  be  t  =  1  for  one  blocklength  CCs.  Then, 

(i)  v,-,  =  (a,  +  a2)~\zu  +  z2l)  iff  (eu,e2i)  =  (0,0). 

(ii)  If  no  error  occurs  at  the  ilh  frame  time,  i.e.,  (eu,e2 «)  =  (0,0),  then 


Vj-i  =  (fli  +  a  2)  \z  u  +  z2i ) 

V  =(fli+  a2)_1[(l+a2)zli  +  (l+al)z2il  +  (al  +  ai)  *[zl*  +  Z2|] 


TIT 


(i)  For  part  (i),  if  (eu,e2 ,)  =  (0,0),  then  by  (23)  and  (24),  one  obtains  v(-1  =  (z1(  +  z2i)(al  +  a2)~‘ 
and 

v,  =  (a  i  +  a2)~x(a2zu  +  ai22 .) 

=  (a,  +  a2r1l(l+a2)zii  +  zu  +(1+0  022.  +  zt.) 

=  (<*i  +  a2)~x[(l+a2)zu  +  (l+ai)2i  ,J  +  (<*i  +  a2yx[zu  +  z2l] 

If  v,'_!  =  (o,  +  a2)_1(zi<  +  z2i)  then  by  (25)  eu  +  e2l  =  0.  Now  since  t  =  1,  and  e2l  cannot 
be  simultaneously  both  nonzero.  In  other  words,  («ii,e2l)  must  be  one  of  the  following  forms: 
(0,0),  (0,o ),  (a,0)  where  a  e  GF( 2*).  Therefore,  it  can  be  readily  verified  that  (eu,e2i)  =  (0,0). 

(ii)  If  {eu,e2i)  =  (0,0)  then  from  part  (i),  it  is  known  that  v,_1  =  (a,  +  a2)~*(z i,  +  z2l)  and  hence  the 
result. 

Note  that  from  (i)  of  Result  1,  it  follows  that  if  (eu,e2i)  =  (0,0)  then 

v,  =  v;-i  +  (*i  +  <z2)”1l(l+a2)zii  +  0+a  l)z2i  1 
Based  on  Result  1,  a  flow  chart  for  the  decoding  algorithm  is  shown  in  Figure  1. 

Example  1: 

A  specific  rate  1/2  dual-3  CC  with  G{D)  =  (1 +D  ,1-KxD )  is  given,  i.e.  a\  =  1,  a2  =  a. 

Solving  v,  and  v,_!  in  Eqs.  (21)  and  (22),  one  has 

v.-i  =  (z  t.  +  z2.)«2  +  («ii  +  el i)°-2 

4 

vi  =  zt.  +  (zi,  +  z2.  a2  +  (e,,  +  e^a2  +  eu 

k 

Then  from  Result  1: 

v,_i  =  (z  i,  +  z  2i  )«2  iff  (e  u  +e2i)  =  (0,0) 


■  v-- 


12 


I -INDEX  OF  THE  RECEIVED 
AND  INFORMATION  SEQ.: 

hi'  hi 

X  -  TRUE:  INDICATING  PREVIOUSLY 
CORRECTED  FRAME 

X  -  FALSE:  INDICATING  PREVIOUS 
FRAME  INCORRECT 

A  -  PREVIOUSLY  CORRECTED 
INFORMATION  SYMBOL 

S  -  TEMPORARY  VARBLE 


(  START  J 


1-  -1 

X-TRUE 

A-  0 

I  -  l  ♦  I 
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and  if  no  error  occurs  at  the  i,k  frame  time,  then 

=  (*ii  +Z2,)ct2 

v.  =zi«  +  (2i,  +22i)«2 

The  detailed  steps  of  the  algorithm  shown  by  the  flow  chart  now  are  described  and  explained  as 
follows: 

Assume  initially  that  the  message  symbol  is  0. 

Step  0  Set  i  =  -1,  x  «  true  and  A  =  0  (the  initial  message  symbol  0). 

Step  1  Set  i  -  j+1,  compute  S  =  (z  tl  +  z^  )a2. 

Step  2:  If  the  previous  time  had  an  error,  i.e.,  x  -  false  then  since  only  one  error  is  allowed 
for  two  continuous  time  frames,  the  current  time  frame  must  be  error  free.  Therefore, 
by  Result  1,  part  (ii),  get  v,_[  =  S  =  (z,,  +  z2l)a2  and  get  v,  =  (zu  +  S).  Next  set 
=  true  and  A  =  v,.  Go  to  Step  1. 

If  the  previous  time  frame  had  no  error,  i.e.,  x  =  true,  then  go  to  Step  3  to  check  if 
any  error  occurs  at  the  current  time  frame  i  or  not. 

Step  3:  If  S  =  A  then  by  Result  1,  part  (i),  leu,e2t]  -  [0,0], 

Thus  v,  can  be  computed  by  Eq.  (24)  v,  =  zu  +  S;  set  A  =  v,  and  x  =  time.  Go  to 
Step  1. 

If  S  *  A  again  by  Result  1,  part  (i),  an  error  occurs  at  time  frame  i ;  set  x  =  false.  Go 
to  Step  1. 

To  illustrate  the  decoding  algorithm  works,  let  the  D-transform  of  input  sequence  be 
x(D )  =  1+c xD  .  Then 

2(D)  =  x(D)G(D)  =  [1,1]  +  [a5,0]D  +  [a,a2]D2 


JlfJt  Jt»A‘Ji'J»' *i'  *•*?»•  ***  tJ*  i^.M  i 
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much  simpler  than  Viterbi  decoding.  In  fact,  less  than  30  gates  are  needed  to  implement  this  decoder. 


VI.  ALGEBRAIC  SYNDROME  DECODING  OF  a  RATE  1/3,  DUAL-K,  CC 


Let  G(D)  =  (l+a]D,l+a2D,l+a3D)  where  a  6  GF( 2*)  is  a  root  of  p(x)  as  defined  in  Section 


V.  Again  in  Eq.  (15),  let 


v(D)=  £  v,D' 
i*0 


2(D)=  X  [zu>Z2i.Z3i]D‘ 

i« 0 


*(£>)  =  z  1,^2. 

i*0 


Substituting  these  equations  and  G(D)  into  Eq.  (15)  and  equating  the  coefficients  yield, 


aiv,-i  +  v,  =  z1(  +  etl 


a2v,-l  +  v.  =  *2,  +  e2i 


a3vi-l  +  v.  =  *3i  +  *3, 


The  solutions  for  v(_,  and  v,  are  given  as  follows: 


For  v,_]  one  has  the  two  solutions. 


vi-\  =  (ZU  +  z3i  +  eU  +  e  3.  Xa  l  +  a3>  1 

v,-i  =  (2  I,  +  22.  +  eU  +  «a)(al  +  a2)_1 


For  v,  one  also  has  the  solutions 


v,  =  aiv.-i  +  +  eu 


v.  =  a2v.-l  +  22,  +  <?2, 


where  v(_,  is  given  in  either  (33a)  or  (33b).  By  (33a)  and  (33b),  one  obtains 


[(a,  +  a2)  1  +  (a,  +  a3r1](2|,  +  *u)  +  (ai  +  a2)  \z2i  +  e2 «) 
=  (a,  +  a3r,(23l  +  e3l) 


Result  2: 


Let  a  rate  1/3,  dual-3  CC  be  generated  by  G(D)  =  (l+a\D ,l+a2D ,1+a^D).  Assume  that  no 
more  than  r  =  2  errors  tan  occur  in  one  constraint  length,  WH{eu>e2i,e  3l)  + 
W///(ei,J+i>c2,i+i>e3,<+i)  ^  2,  where  denotes  Hamming  weight 

(i)  If  [e  lf  ,e  2 i  ,e  3l  ]  =  [0,0,0],  then 

[(a,  +  a2)~l  +  (a,  +  aYTX]zu  +  (*i  +  a2)~xz2i  =  (aj  +  a3r‘z3l  (36) 

vt_i  =  (z 1,  +  z  2l  )(a  1  +  a2)_1  (37) 

V;  =  V;_!  +  Z  j,  (38) 

(ii)  If  [(a,  +  a2)~x  +  (a,  +  a3)-1]zll  +  (a j  +  a2)~xz2i  =  (a,  +  a3)~'z3j ,  then  at  least  one 

error  occurs  at  frame  i,  i.e.,  3l)  S  1. 

(iii)  If  [(a,  +  a2)_1  +  («i  +  a3)_1]zu  +  (dj  +  a^z^  »  {ax  +  a3)_1z3l,  then  either  0  or  2 
errors  occur  at  frame  i,  i.e.,  W/w(eli,e2l,e3l)  =  0  or  2. 

(iv) 

v.-l  =  (zl,  +  z2«)(°l  +  a2)~l  (39) 

and 

[(a,  +  a2yl  +  (a,  +  a3)_1]zi,  +  (ai  +  02)  *z2i  =  (al  +  al)  *z3i  (40) 

iff  [e1,,e2l,C3,]  =  [0,0,0], 

(v)  Suppose  [(a,  +  a2)-1  +  (aj  +  a3)-1]z1(-  +  (a]  +  a2)-1z2l  *  (a,  +  a3r'z3l  and 

WH(eu<e2i>ex)  =  1- 


(a)  If 


a !  +  a2 
07  +  a3 


z2i 


di  +  flj 

-  z 3i (a  1  +  fl2)v.-l  =  0, 

a2  +  a3 


(41) 


Qi  +  <*2 
a2  +  a3 


then  e1(  =  zu  + 


a,  +  a3 


z2i  + 


z3i  *  0,  e2l  =  e3(  =  0. 


CL  i  +  Cly  CL  i  +  (2  2 

(b)  If - - -  z2i  + - - -  z3i  +  (ax  +  a2)v,_ ,  *  0  then  eu  =  0. 

CL 2 +  ^2  *  ^3 

Proof: 

(i)  Follows  immediately  from  Eqs.  (33),  (34)  and  (35). 

(ii)  Follows  immediately  from  part  (i). 

(iii)  If  [(fl!  +  a2)_1  +  (a,  +  a3)'!]zu  +  (al  +  a2)~lz2i  =  (a1  +  a3)'z3(,  then  by  (35)  one 
obtains 

[(a,  +  d2)-1  +  (a  i  +  a3r']«ii  +  («i  +  =  ( a\  +  a3)”'«3i  (42) 

Now  if  (e  i,  ,e ^  ,e  3l )  =  1,  then  it  can  be  readily  verified  that  (42)  cannot  be  satisfied.  Moreover, 
since  for  dual-k  of  rate  1/3  CC,  the  correcting  capacity  t  =  2,  one  concludes  that  Ww(ell,e2l,e3l)  =  0 
or  2. 

(iv)  If  v;_!  =  (zu  +  z2i)(al  +  a2)~l  and  {(oj  +  a2)~'  +  (a,  +  a3r']zi,  +  (fl]  +  a2)~'z 2i  = 
(a,  +  a3)_1z3(,  then  by  (33),  (35),  (38)  and  (40),  one  obtains 

eu  +  =0  (43) 

[(at  +  a2)_1  +  (a,  +  a3r']eu  +  {ax  +  a2)~xe2i  =  (a,  +  a3)_Ie3l  (44) 

A  substitution  of  (43)  into  (44)  gives  eu  =  e3i.  This  implies  eu  =  e %  =  e3l.  Since  t  =  2,  three 
errors  are  not  allowed.  Thus,  (eii,e2i,e3i)  =  (0,0,0).  The  converse  follows  immediately  from  part  (i). 

(v) 

(a):  A  substitution  of  (41)  into  (33b)  yields 

(a,  +  a  2)(a  2  +  a3r‘(z2l  +  z3()  =  z ,,  +  z2,  +  «i,  +  (45) 


Solving  (45)  and  (35),  one  obtains 
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eh  =  «3i 

a  \  +  q-$  + 

el<  =  zh  +  - — : -  *2i  +  ~ — ~~~  z3i  +  «3. 


+  a-i 


a 2  +  *3 


(46) 

(47) 


Since  VVH(ej,  =  1,  e2l  and  e 3l  cannot  both  be  nonzero  thus,  ell=e-i=0  and 


fli  +  ai 


el<  ~  zl<  + 


^2  ■*"  ^3 


z2i  + 


a !  +  a2 
a2  +  <23 


Z-Ki  *  0. 


(v) 


(b):  Eliminating  v,  in  Eqs.  (30)  and  (31)  gives 


eU  +  e2 1  +  zl,  +  22.  +  <*1  +  fl2vi-l  =  0. 


(48) 


Similarly,  eliminating  v,  in  Eqs.  (31)  and  (32)  gives 

eU  +  e3 .  +  zl.  +  z3i  +  al  +  a3vi-l  =  0-  (49) 

Adding  Eqs.  (48)  and  (49)  yields 

e2i  +  3ei  +  z2i  +  z3i  +  o2  +  a3v,_,  =  0 
or 


a  1  +  <2  2  j  +  fl  2 

777  («2<-  +{i)  =  77-  (*2i  +z3i)  +  «i  +«2Vi-i  *  o. 

Thus  e2<  *  e 3i •  Since  VV#  (ej,,^,,^)  =  1,  one  obtains  e1(  =  0. 

The  flow  chart  for  decoding  a  rate  1/3,  dual-3  CC  is  shown  in  Fig.  2.  Initially,  i  =  -1  and 
A  =  0,  where  i  is  the  index  for  both  received  and  information  sequence  and  A  is  used  to  store  the  pre¬ 
viously  corrected  message  symbol.  For  any  two  contiguous  time  frames  (time  frame  i  and  i  +  1),  four 
Boolean  variables  BOLN  [j  ],  1  <  j  <  4  are  defined  as  follows: 

If  ((a]  +  a2)_l  +  (ai  +  <23)_*  ]z ii  +  (a i  +  fl2)_* lz 2i  =  (^l  +  03)23,  then  BOLN[  1J  =  true, 


else  BOLA/fl)  -  false. 


STAIT 
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If  [(a,  +  a2)~x  +  (fli  +  a3r1]zlti+1  +  fat  +  a2)_1]z2.i+1  =  (a\  +  a3)2  3.1+i  then  BOLN\2\ 
«  true,  else  BOLN  [2]  =  false. 

If  (zj,-  +  z2t)(al  +  a2)_1  =  A  ,  then  BOLN[3]  =  true,  else  BOLN[2]  =  false. 

If  — — —  r2i  +  — - - -  z3i  =  (aj  +  a2)v1_1  then  BOLN[4)  =  true,  else  BOL\\4\ 


<*2  +  a3 

false. 


a  2  +  a3 


Four  cases  are  discussed  according  to  the  Boolean  values  of  BOLN  [1]  and  BOLN  [2], 

Case  1:  B0LN[1]  and  BOLN  [2]  are  both  true. 

By  Result  2,  part  (iii),  there  are  0  or  two  errors  for  both  time  frames  i  and  i  +  1,  i.e. 
Wti(eu,e 2i,e3i)  =  0  or  2,  WH(elii+1,«w+„e3fi+1)  =  0  or  2. 


If  BOLN[3]  is  true,  then  by  Result  2,  part  (iv),  (e  1i,«2<,e3i)  =  (0-0,0),  and  thus  by  Result  2,  part 
(i),  v,  can  be  determined  by  Eq.  (38). 

If  BOLN[ 31  is  false,  then  by  Result  2,  pan  (iv),  VV'w(e]l,e2(,e3l)  >  1,  and  since 
Wlf(e  i,  ,eh  ,e3l )  =  0  or  2,  as  discussed,  one  concludes  that  WH(eu  ,e2l  ,e3()  =  2.  Thus 
WH(e ,  1+1,e2l+1,e3  l+1)  =  (0,0,0)  for  the  reason  that  no  more  than  2  errors  can  occur  in  two  contiguous 
Dtne  frames.  Therefore  v,  can  be  determined  by  Eq.  (37)  v,  =  (z1(+1  +  Z2,i+i)(ai  +  a2)  • 

Case  2:  BOLN  [  1 1  false,  BOLN  [2]  true. 

By  Result  2,  pan  (ii),  WH(e  1(  ,e2l  ,e3l )  >  1.  By  Result  2,  pan  (iii),  VVw(e ,  1+1,<?2  )+1,e3 1  +  1)  =  0  or 
2.  Also  by  the  assumption  that  WH(e] ,)+1,e2,;+I,e3tl+|)  +  Ww(e i,»+i-e2,i+i-e3.«+i)  -  2.  Thus, 
Wn(t\.,  +  \’*2J*\'*\,*\)  =  0'  *e.  (e  i,i+i 2,i +t,e 3,i+t'  z  (0-0-0)-  and  so  v,  again  can  be  determined  by 
Eq.  (37). 


Case  3:  BOLN\\  \  true,  BOLN [2]  false. 
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By  Result  2,  part  (iii),  W H  (e  u  >e  %  ,e  n)  =  0  or  2.  By  Result  2,  part  (ii), 
wh(* i.«+i.c2^+i.«3,i+i)  -  Also  by  the  assumption  that  («,,,<? 2.,* 3.)  +  Ww(e,  1+l,e2(+1,e3(+l)  -  -• 

Thus,  WH(eii,e2i,eu)  =  0  and  so  by  Result  2,  part  (i),  v,  can  be  computed  by  Eq.  (38). 

Case  4:  BOLN[\]  false,  BOLN[ 2]  false. 

By  Result  2,  part  (ii),  W’tf(e,i,e2i,e3i)  £  1  and  M/w(eu+1,e2ii+,,e3j+1)  <  1  and  since 
w'//(«u.«2i»tf3i)  +  w'tf(«i,.+i.«w+i.«v+i)^  2,  it  follows  that  WH(eu,e2i,e3i)  =  l  and 
^//(ci,i+i>e2,i+i*e3,i+i)  =  1-  Therefore  by  Result  2,  part  (v),  if  BOLN[4]  is  true  then  e2l  =  0,  other¬ 
wise  eu  =  0.  And  so  v,  can  be  computed  by  Eq.  (34b)  or  (34a)  depending  on  e2l  =  0  or  eu  =  0. 

After  v,  is  determined,  then  v,  is  substituted  for  A  and  the  recursive  process  is  repeated  to  deter¬ 
mine  v,-+|.  Finally  the  entire  estimated  information  sequence  is  obtained.  The  estimated  information 
sequence  is  exactly  the  original  input  information  sequence  as  long  as  there  is  no  more  than  2  errors 
for  each  constraint  length. 

To  illustrate  the  above  algorithm,  an  example  is  given  as  follows: 

Example  2: 

Consider  a  special  1/2  rate  dual-3  CC  with  G(D)  =  (l+D,l-KxD,l+a6£>),  i.e.,  A,  =  1,  A2  =  a, 
A  3  =  a6  then  the  four  Boolean  functions  are  as  follows: 

Bolu{  1]  =  (a4z,,  +  a6z2,  =  z  3| ) 

Bolu[  2]  =  (a4zli+1  +  a6zi)+1  =  Z3)l+i) 

Bolu[  3]  =  ((zu  +  z^Ja2  =  A ) 

Bolu{  4]  =  (a  3z,,  +  a3z3(  =  v,_,a5) 

Let  the  input  be  x(D)  =  1  +  aD2.  Then  the  encoder  output  is 

y(£>)  =  x{D  )G(D)  =  [1,1,1]  +  [l,a,a6]£>  +  [a,a,a]D2  +  [a,a2,l]D3 


U  *  V  'T.*i  V 


Lii 
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Next  let  e_{D)  =  [a  ,0,0]  +  [0,0, a]D  +  [a,ot2,0]D3,  so  that  the  received  sequence  z(D)  =  (a, 1.1) 
+  [  1  ,a,a2]D  +  [a,a,a]D2  +  [0,0, 1]D3. 

The  decoding  algorithm  then  can  be  checked  as  follows: 

Step  0:  i  =  -1,  A  =0. 

Step  1 :  i  =  i+1  =  0. 

Step  2:  a4z10  +  a6z20  =  a3  *  1  =  230  BOLN[  1]  is  false. 
a4zu  +  a6z21  =  a6  *  a2  =  z31  ->  BOLN[ 2]  is  false. 
a3z20  +  a3Z3o  =  0  =  v_j  -4  BOLN[4]  is  true. 

Step  3:  el0  =  z10  +  a2z20  +  a3zM  -  Ot5,  v0  =  z  10  +  v_,  +  el0  «  a  +  0  +  a5  =  1  and 
A  =  v0  =  1. 

Step  1  i  =  i+1  =  1. 

Step  2:  a4z„  +  a6zZ(  =  a6  *  a2  =  z3,  -+  BOLN[  1]  is  false. 
a4z  12  +  c^z^  =  a  =  z32  -4  BOLN [2J  is  true. 

Step  3:  v,  =  (zl>1+1  +  z2<1+,)Ct2  -  (z1<2  +  zu)a2  =  0,  A  +  v,  =  0. 

Step  I  i  =  i  +  1  =  2. 

Step  2:  ct  Z|2  +  ct 6z22  =  oc  =  z32  -4  BOLN  [1]  is  true. 

a4z13  +  a6z23  =  0  *  1  =  z33  -4  BOLN\ 2]  is  false. 

Step  3  v2  =  v ,  +  z  12  =  0  +  a  =  a 

Step  I :  i  =  i+1  =  3. 

Step  2:  a4z13  +  a6z23  =  0  *  1  =  z33  -»  BOLN[  1)  is  false. 
a4Z]4  +  a6z24  =  0  =  -4  BOLN [2]  is  true. 


Step  3:  v3  =  (z ,  4  +  z2>4)CC2  =  0,  A  =  v3  =  0. 

Step  1  i  =  i+1  =  4. 

Step  2:  a4z14  +  a6z24  =  0  =  z^  — >  BOLN[  1]  is  true. 
a4z15  +  a 6zm  =  0  =  z35  — *  BOLN[ 2]  is  true. 

(z  14  +  z24)a2  =  0  =  A  — *  BOLN[3]  is  true. 

Step  3:  v4  =  v3  +  z  14  =  0. 

Finally,  v0  =  1,  v ,  =  0,  v2  =  a,  v3  =  0  and  v4  =  0  or  v(D)  =  1  +  a D2,  which  is  the  same  as  the 
original  input  sequence.  A  simulation  program  for  the  algorithm  was  written  in  Pascal.  An  LSI  design 
of  the  algorithm  is  given  in  Appendix  B. 

The  implementation  of  this  decoder  is  much  simpler  than  Viterbi  decoding,  thus,  this  decoder  can 
be  applied  to  telephony  and  HF  radio,  where  a  moderate  error  correction  is  desired  at  a  relatively  low 
cost,  see  [7], 

V.  CONCLUSION 

In  this  report,  algebraic  syndrome  decoders  of  both  rate  1/2  and  1/3,  dual-k,  nonsystematic  con¬ 
volution  codes  are  developed  in  detail,  including  an  efficient  method  for  finding  the  corrected  message 
sequence.  This  is  achieved  without  the  necessity  of  storing  a  number  of  states  or  paths  in  a  constraint 
length  of  the  error  trellis  diagram.  Currently,  the  problem  of  algebraic  syndrome  decoding  of  a  rate 
1  In.  dual-k,  CCs  for  n  >  3  is  being  investigated.  Finally,  a  single  VLSI  chip  with  CMOS  technology 


for  the  rate  1/2,  dual-2  algebraic  syndrome  decoder  is  being  developed. 
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